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ABSTRACT 

Axisymmetric magnetohydrodynamic (MHD) simulations have been applied 
to investigate the interrelation of a central stellar magnetosphere and stellar 
wind with a surrounding magnetized disk outflow and how the overall formation 
of a large scale jet is affected. The initial magnetic field distribution applied 
is a superposition of two components - a stellar dipole and a surrounding disk 
magnetic field, in both either parallel or anti-parallel alignment. Correspondingly, 
the mass outflow is launched as stellar wind plus a disk wind. Our simulations 
evolve from an initial state in hydrostatic equilibrium and an initially force-free 
magnetic field configuration. Due to initial differential rotation and induction of a 
strong toroidal magnetic field the stellar dipolar field inflates and is disrupted on 
large scale. Stellar and disk wind may evolve in a pair of collimated outflows. The 
existence of a reasonably strong disk wind component is essential for collimation. 
A disk jet as known from previous numerical studies will become de-collimated 
by the stellar wind. In some simulations we observe the generation of strong 
flares triggering a sudden change in the outflow mass loss rate by a factor of two 
and also a re-distribution in the radial profile of momentum flux and jet velocity 
across the jet. We discuss the hypothesis that these flares may trigger internal 
shocks in the asymptotic jets which are observed as knots. 

Subject headings: accretion, accretion disks - MHD - ISM: jets and outflows - 
stars: mass loss - stars: pre-main sequence galaxies: jets 

1. Introduction 

Astrophysical jets are launched by magnetohydrodynamic (MHD) processes in the close 
vicinity of the central object - an accretion disk surrounding a protostar or a compact 
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objec t (IBlandf ord fe Pavnd Il982l ; iPudritz fe Norman! Il983l ; ICamenzindl Il99dl ; iPudritz et al. 



20071 ; ICabritll2007l ). However, the principal mechanism which actually launches the outflow 
- the transition from accretion to ejection - for a certain disk at a certain time is still not 
completely understood. 

During the recent decade, numerical simulations of MHD jet formation became more 
and more feasible and substantially helped to improve our understanding of how jets emerge. 
In general, these simulations may be disti nguished in tho s e taki n g into account t h e tim e 
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The first approach allows to study the launching process directly in particular the mech- 
anism lifting matter from the disk plane into the outflow. However, this approach is compu- 
tationaly very expensive and yet limited by spatial and time resolution. In order to study the 
jet formation process - the acceleration and collimation of a disk/stellar wind - it is essential 
to follow the dynamical evolution of the system for many (several thousands) of rotational 
periods and on a sufficiently large grid with appropriate resolution. For such a goal, the 
second approach is better suited. Since it is computationally less expensive it also allows 
to perform series of simulations for parameter studies. It is clear that the prescription of 
mass flow rate and magnetic flux profile constrains the result of the simulation more than a 
consistent simulation of the jet-disk evolution which could in fact provide the mass loss rate 
from the disk into the jet. On the other hand, the current status of MHD disk modeling has 
its own limitations. In particular the magnetic field structure in the disk is a rather open 
question unless radiative MHD, global simulations of dynamo-active disk models provide 
fully self-consistent results. The aim of this paper was to investigate the interaction and 
interrelation between a stellar wind and a disk wind and represents a unique approach in 
that field. For this first step it is of advantage to govern the simulation by well understood 
boundary conditions. Further, for our goal it is essential to study the long term evolution of 
the outflow at considerable distances from the star. If we would include the disk evolution 
in the simulation it would be hard to reach the appropriate time scales. In fact all disk 
simulations so far stop at earlier time scales. Future work should include the disk evolution 
for the jet launching. 

One may further distinguish between the different i nitial setup for the s e simulation s 
- some of them consider a pure stellar dipole (see e.g. lUchida fe Shibatal (119841 . Il985l ); 
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opc 

f l2003f ): lFendtl (120061 )). The case of a superposed stellar and disk magnetic field is yet rarely 
tre ated in simulat i ons, h owever, the first such model configura tion was discussed already 
by lUchida fe Low! (I198ll ) . One example is iMiller fe Stond (119971 ) who superposed a central 
dipole with an aligned vertical disk field. 

In this paper, we study the long-term evolution of a two-component MHD outflow 
consisting of a stellar wind launched from a stellar magnetosphere and a surrounding disk 
wind. It essential to follow the time evolution for very long term in order to be able to take 
into account the evolution also of the outer regions of the disk magnetosphere as much as 
possible. 



Goodson et all (119971): Miller fe Stonei (119971) : lFendt fe Elstner! fll999l.l2000h:lRomanova et al 



d2002h: iMatt fe Pudritzl (hood)), others a pur e disk field (see e.g. 
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2. Stellar magnetospheres inside disk outflows 

The magnetic field of protostellar jets (and probably also those of micro-quasars) most 
probably consists of two components - a central stellar, probably dipolar field plus a field 
component provided by the surrounding accretion disk (either generated by a disk dynamo 
or advected from interstellar space or both). 



2.1. Impact of a stellar magnetosphere on the large-scale outflow 

In the following we qualitatively expose the main aspects of the interaction between 
a central stellar field and the disk magnetic field and how that may affect the overall jet 
formation. 

Enhanced magnetic flux. The stellar field adds magnetic flux to the system. Assuming 
a polar field strength B* and a stellar radius R±, the dipolar field scales with 

This has to be compared to the disk poloidal magnetic field provided either by a disk dynamo 
or by advection of ambient interstellar field. The latter can be estimated by equipartition 
arguments, and is limited to 
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The stellar magnetic field will not remain clo s ed, bu t will inflate due to shear between the disk 



surfac e and the star (e.g . [Uchida fc Shibatal (119841 ) ; iFendt fc Elstnerl (120001 ) ; lUzdensky et al. 



(120021 ); iMatt fc Pudritzl (120051 ) ) . The additional Poynting flux that threads the disk may 
support the MHD jet launching and may provide an additional energy reservoir for the 
conversion of magnetic energ y to jet kineti c energy, thus, implying a g re ater asymptotic 
j et sp eed (Michel scaling, see iMichell (119691 ) : IFendt &: Camenzindl (119961 ); IFendt fc Ouyed 
J2004K 



Additional magnetic pressure. The central magnetic field provides additional magnetic 
pressure implicating possible de-collimation of the overall outflow. The stellar magnetic field 
may drive a strong stellar wind which will remove stellar angular momentum. This stellar 
outflow will interact with the surrounding disk wind. The observed protostellar jets may 
consist of two components - a stellar wind and the a wind, both with a strength depending 
on intrinsic (yet unknown) parameters. 

Angular momentum exchange. In the scenario of magnetic "disk locking", the stel- 
lar field threading the disk will re- arrange the global angular momentum budget. Angular 
momentum of the star is transfered by the dipolar field and deposited at the inner disk. 
Thus, matter orbiting in this region will be accelerated to slightly super-Keplerian rota- 
tion. This has two interesting implications, (i) Due to the super-Keplerian speed this disk 
material could be more easily expelled into the corona by magneto-centrifugal launching 
(IBlandford fc Paynd Il982l ; iFerreiral 119971 ) and forms a disk wind, (ii) Excess angular mo- 
mentum in that disk area will slow down accretion unless removed by some further process. 
The disk outflow launched from this very inner part of the disk can be an efficient way to do 
this. This scenario is similar to the X-wind models (jShu et al.lll994t iFerreira et al.l 120061 ). 



Non-axisymmetric effects. In addition to the simple picture of an axisymmetric config- 
uration, an inclined stellar magnetic dipole will add non-axisymmetric effects. A moderate 
non-axisymmetric perturbation may result in warping of the inner disk, and, thus, a preces- 
sion of the outflow launched from this area. For extreme cases of inclination, jet formation 
may be completely prevented. 

A rotating inclined dipole further implies a time- variation of the magnetic field strength 
at the inner disk radius. This may lead to a time- variation in the accretion rate and, also, 
the m ass outflow flow rate. Numerical simulations of the warping process (jPfeiffer fc Lai 
20041 ) indicate that the warp could evolve into a steady state precessing rigidly. Disks may 
be warped by the magnetic torque that arises from the a slight misalignment between the 
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disk and star's rotation axis 



2.2. Simulations of star-disk magnetospheres 



The simulation of outflows from star-disk magnetospheres is a difficult numerical task 
due to strong gradients in magnetic field strength, density/pressure, and also the velocity 
field. Still, the first simulations of a stella r dipole surrounde d by a disk have been pre- 
sented in seminal papers as early as 19 84 bvlUchida fc Shibatal (1 1984. 1985 ). Probably due 
to the success of the disk wind models ( Blandford fc Payne! 1982 ; Pudritz fc Norman 1983 ) 
the formation of outflows from stellar magnetospheres became somewhat unatte nded before 
the topic was re-disco vered in the early 90ies with several model suggestions (jCamenzind 
199CM : IShu et al.lll994j ). Numerical, stationary- state MHD solutions of the star-disk jet for- 



mation problem were pre sented subsequently (IFendt et al.l Il995l ; iFendt fc Camenzindl Il996 
Sauty fc Tsinganoslll994l ). 



Ti me- dep endent simul a tions of dipolar magnet ospheres started again later lHayashi et al. 



(1996 


); 


Hirose et al. 


(1997 


); 



lations including the disk evolution was short - a few inner disk rotations. 



Goodson et al. (IGoodson et al.l 119971 . Il999l ; iMatt et al.l |2002| ) in a series of papers 
presented simulations runs of dipolar star-disk magnetospheres for up to 150 (inner) disk 
rotations, in particular considering the reconnection/flaring behavior of the magnetic field 
close to the inner disk radius. As another result, a highly collimated, narrow jet emerged 
along the rotat ion axis for which the authors also presented forbidden emission line maps. 
Simulations by IMatt et al.l (120021 ) followed a similar setup, however, based on an increased 
grid scale of (0.8 x 3.2AU) and run for times scales of 150 days (stellar case). Magnetic flares 
emerging from the inner disk location lead to episodic mass ejections. These flares were 
both triggered and triggering a variation of the disk mass accretion rate and a corresponding 
change in the inner disk radius. Reconnection is strongly interrelated to magnetic diffusivity 
(resp. electric r esistivity), howeve r, such model assumptions were not particularly specified 
in these papers. IMatt et al.l (120021 ) investigated the setup of an ambient vertical field aligned 
with a central dipolar field. Their scenario is somewhat different to what we consider: the 
surrounding magnetic field is rather weak and vertical and, thus, not an active outflow- 
generating disk field as in our case. The authors identified three classes of field geometries 



1 This disk warping mechanism may also operate in the absence of a stellar magnetosphere as purely 
induced by the interaction between a large-scale magnetic field and the disk electric current and, thus, may 
lead to the precession of magnetic jets/outflows (Lai 20031 ) 
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Fig. I. — Scheme of model setup. Along the equatorial plane (r-axis) we distinguish jet 
inflow boundary conditions (shaded areas of ghost cells outside active grid) from the star 
(r = 0.0, r* = 0.5) and the disk region from r« = 1.0, r = r out ). The mass flux from star 
or disk is prescribed by the profile for the inflow density and inflow velocity (dashed arrows). 
The latter is typically of about 0.1% of the local Keplerian speed. The stellar boundary is 
in rigid rotation, the disk boundary in Keplerian rotation. The disk inner radius is at the 
co-rotation radius. In the gap area between disk and star (r = 0.5, ...,Tj = 1.0) a minimum 
"floor" mass flux is defined in order to keep the simulation going. Thin contours indicate 
initial poloidal field lines for one of the simulation runs. 
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- as in our setup depending on the direction of the ambient field. As a main result, in all 
th ree cases the collimation of the central outflow was observed, similar to earlier suggestions 
by lKwan fc Tademarul (119881 ). 



In comparison, iFendt fc Elstnerl (119991 . 120001 ) did not treat the disk structure in their 
simulations, allowing them to follow the evolution of the the stellar magnetosphere for more 
than 2000 rotation periods of the inner disk. These simulations sho wed for the first time 
that the axial jet feature o bserved in some of the simulations before (IMiller fc Stone! 119971 ; 
Goodson et al.l 119971 . Il999l ) is in fact an intermittent feature and disappears on the long 
term. The initial axial jet appears essentially during rel axation process of the hydrostatic 
initial condition towards a hydrodynamic steady state. IFendt fc Elstnerl (120001 ) were first 
to introduce a physical stellar wind boundary condition in addition to the disk wind. The 
emerging two component outflow from star and disk remained un-collimated. This was 
understood as caused by the fact that no net electric current was driven along the outflow 
(since the initial dipolar field distribution). Since the stellar dipolar field becomes very weak 
along the disk, a classical magneto-centrifugal wind launching becomes inefficient. This is 
another reason for the lack of collimation observed in these simulations. 



Simulations of the disk-star interaction iRomanova et al.l (120021 ); iKuker et al.l (120031 ); 
Romanova et al.l (120041 ) focus on the dipolar accretion process and potential angular mo- 
mentum exchange between star and disk. No evidence was found for persistent outflows 
and only episodic ejections from the inner disk area occurred. Time-dependent simulations 
of winds from an initially sp herically symmetric stellar magnetic field were performed by 
Keppens fc Goe dbloed (1999), evolving exactly into the classical stationary Weber-Davis 



solution. iMatt fc Pudritzl (120081 ) derive the angular momentum loss of pure stellar winds 
by numerica l simu lation s. The setup is similar to th e stationary-state models calculated by 
Fendt et al.l (119951 ) and IFendt fc Camenzindl (119961 ) . The outflow angular momentum loss 
derived from these simulations is comparable to what can be derived from the stationary 
models^. 

While the simulations from dipolar magnetospheres failed to show collimated outflows, 
MHD simulations of disk winds did actually proof the self -colli mating characteristics o f 
MHD winds. After seminal papers by lUstyugova et al.l (119951 ) and lOuyed fc Pudrita (119971 ) 
this approach was fu rther developed taking in to account the time-dependent change of 
disk field inclination (|Krasnopolsky et al.l 119991 ). turbulent magnetic diffusivity in the jet 
( IFendt fc Cemeljid 120021 ). a variation in the disk wind magnetic field and density profile 



2 Sean Matt , priva t e communication. Number values for the angular momentum flux were not included 
in lFendt et all |l995f ): IFendt fc Camemrindl |l996h 
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dPudritz et al.lbood: [j^dtll2006h . or non-axisymmetric effects (jOuyed et all2003l;lKigure Sz Shibatal 



2005 



Anderson et al 



Casse fc Keppens 



2002 



2006 ), and also the disk dynamical evolu tion (jKudoh et al.l 



f998 



Kigure fc Shibatal 120051 : IZanni et al.l 120071 ). 



Recent simulations by iMatsakos et al.l (120081 ) investigated the "topological stability" of 
two-component (star-disk) self-similar solutions derived from stationary MHD. In difference 
to our approach their simulations do not start from an initial hydro-sta^c state, but from an 
initial dynamical steady state solution of the MHD equations, with analytical extrapolations 



i n the case of radially self-similar solutions (which are singular on the axis). iMeliani et al. 



( 120071 ) presented two-component outflow simulations including the treatment of the disk 
evolution. However, although the wind dynamics indeed consists of two components, the 
magnetic field distribution basically resembled a monotonous field profile (i.e. no stellar 
dipole, no stellar outflow involved). 



3. Simulation model setup 



We perform axisymmetric MHD simulations of jet formation for a se t of different mag- 
netic fi eld geometries and mass fluxe s. The gener al model setup follows lOuyed fc Pudritz 
( 119971 ) . iFendt &: Cemeljid (120021 ) and iFend t (2006 ), however, with i mportant modifications . 
The original ZEUS-3D ideal MHD code Jstone fc Normanl [l992a1 Jbl: lHawlev fc Stonlll995l) 
extend ed for physical magnetic resistivity (see description and tests in IFendt fc Cemeljic 
( 120021 ) ) is used. For the purpose of this paper the magnetic diffusivity was set to such a 
low level that it does not affect the overall collimation of the outflow. However, resistiv- 
ity/ diffusivity is essential for our simulations, as it allows for magnetospheric reconnection 
phenomena. 



The set of MHD equations considered is the following, 



dp 
di 

du 



+ V • (pv) = 0, V • B = 0, —j =VxB, 



+ (v ■ V) v 



+ V(p + p A ) + pV$ - 



jxB 



— - V x [vx B r)j 

at \ c 



P 



An 



+ p(V ■ V) jT)j 



0. 



(3) 
(4) 
(5) 
(6) 



with the usual notation for the variables (see IFendt fc Elstnerl (120001 ); iFendtl (120061 )). The 
magnetic diffusivity is space and time dependent and is denoted by i](r, z). 
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We a pply a polytropic equa tion of state for the gas with the polytropic index 7 = 
5/3. As (lOuyed fc Pudritzl 119971 ) we have added turbulent Alfvenic pressure in order to 
allow to keep the corona "cool". In fact, we do not solve the energy equation [H] and apply 
the internal energy of the gas reduced to e = p/(j — 1). Two major reasons to do this 
are both computational speed and numerical stability. Some of our long-term simulations 
did last already more than two months of CPU time and it would be impossible to reach 
the desired evolu t ionary steps if the energy equation would have been solved. Following 
Ouyed fc Pudrita (119971 ) this approach also allows to combine gas pressure forces and gravity 
under the same derivative in the code. Thus, instead of subtracting gradients, we apply 
the gradient of the difference, which results in perfect stability of the initial state. Note, 
however, that recent work by Clarke and collaborators seems to indicate that relaxation of 
the polytropy assumption may affect the dynamical evoluti on in certain domains of the jet, 
in particular regions with shocks or contact discontinuities ( IRamsey fc Clarke 1 120041 ) . 



The magnetic diffusivity can be considered as turbulent and, thus, be related to the 
Alfvenic turbulent pressure pa if we assume that it is primarily due to the turbulent Alfvenic 
waves which are res ponsible for the turbulen t Alfvenic pressure applied in the simulations. In 
our previous work (IFendt fc Cemeljidl2002l ) we have derived a toy parameterization relating 
both effects by parameterizing the turbulent magnetic diffusivity similar to the Shakura- 
Sunyaev parameterization of turbulent viscosity, rjt = a m vl, where a m < 1 and I and v are 
the characteristic dynamical length scale and velocity of the system, respectively. With that 
we obtain i] t = a m v t l and with f3 t = (c s /v t ) 2 and = jp/p, it follows that vl = 
or, normalized, v'\ = ^p' 7-1 . For the chosen polytropic index this implies a magnetic 
diffusivity 77 ~ p 1 / 3 if I is constant (see also the discussion in iFendtl ( 20061 )). We again 
stress the point that in the present pa per the magnetic d iffusiv ity is on a such a low level 
(much below the critical level found in lFendt fc Cemeljid (120021 )) that it does not affect the 
dynamics, but allows for reconnection. 

For the numerical grid, we use the "scaled grid" option by ZEUS with the element size 
decreasing inwards by a factor of 0.99. The size of the cylindrical grid is (256 x 256) elements 
resulting in a physical grid size for all simulations is (r x z) — (80 x 80)r; corresponding to 
(4 x 4) AU for r; ~ 10 R Q . Thus, the disk gap (r < 1) is resolved with 11 grid elements. 
Time is measured in rotation periods (Keplerian orbits) at the inner disk radius. 



I n summary, compared to our previous studies (IFendt fc Elstnerll2000l ; IFendt fc Cemeljic 



20021 ; IFendtl 120061 ) . the main new feature included in the present approach is that the initial 
magnetic field distribution consist s of two components - a stellar dipolar field and a disk field. 
In difference to recent studies by (jMeliani et al.ll2007l ) we included both outflow components 
and magnetic field components in the simulation box, in particular treating the star-gap-disk 
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boundary. 

3.1. Boundary conditions 

Along the r-boundary we distinguish between the star extending from r = 0.0 to r = 
= 0.5, the gap region extending from r = 0.5 to r = = 1.0, and the disk region 
from r = 1.0 to r = r out , see FigJTJ A Keplerian disk is taken as a (fixed in time and 
space) boundary condition for the mass inflow from the disk surface into the corona and the 
magnetic flux. The stellar surface is approximated by a rigidly rotating "disk" in cylindrical 
coordinates. The stellar rotation is chosen such that the inner disk radius is located at the 
co-rotation radius. 

The initial magnetic field is purely poloidal. Magnetic field lines are anchored in the disk 
and the rotating star and are in co- rotation with their respective foot point. The poloidal 
magnetic field profile along the r-boundary remains fixed in time and is, hence, determined 
by the choice of the initial magnetic field distribution. 

The disk region governs the mass inflow from the disk surface into corona (denoted 
as "disk wind"). In addition to that we prescribe a stellar wind with different mass load. 
The hydrodynamic boundary conditions are "inflow" along the r-axis for the disk and stellar 
region and either "inflow" (very light mass flow) or "reflecting" for the gap region, "reflecting" 
along the symmetry axis and "outflow" along the outer boundaries. Matter is "injected" 
from the disk and the star into the corona parallel to the poloidal magnetic field lines 
with very low velocity v in j(r, 0) = ViV^{r)B^>/ Bp and with a density p in j(r, 0) = T]ip(r, 0). 
The proportionality constants are typically v x ~ 10~ 3 and r\\ ~ 100 for both, stellar and 
disk wind (but usually different for both components). Along the gap we impose a floor 
value for the density r\\ ~ 10~ 3 and a similarly low velocity. For the injection velocity, the 
assumption is that the initial disk wind speed is in the range of the sound speed in the disk, 
Vinj{r) ~ c s (r) ~ v kcp ~ r" 1 / 2 . 

These mass loss rates from star and disk, respectively, are our other main parameters 
besides the respective magnetic flux (see Tab. [[]). 

3.2. Initial conditions 

As initial state we prescribe a force-free magnetic field along with a gas distribution 
in hydrostatic equilibrium. Both is essential in order to avoid artificial relaxation pro- 
cesses caused by a non-equilibrium initial condition. The initial density distribution is 
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Fig. 2. — Example initial magnetic field distributions. Shown are poloidal magnetic field 
lines. Full and dashed lines indicate the direction of magnetic flux. The plots show different 
options for the strength and orientation of the superposed stellar and disk magnetic field 
components ^o.d = 0.01,-0.01,-0.1, resp. ^ ,* = 5.0,5.0,3.0 (from top to bottom; runs 
A2, A3, A4a). 
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p(r,z,t = 0) = (r 2 + z 2 )~ 3 / 4 . As initial magnetic field distribution we superpose a dipo- 
lar steJJ_ai^^ld_and a disk poten tial field. For the disk fie ld component we apply the model 
of lOuyed &: Pudrita (119971 ) and iFendt fc Cemeljid (120021 ) . We prescribe the magnetic field 
distribution as derivative of the magnetic flux distribution \I/(r, z) = J B p dA (i.e. the 0- 
component of the vector potential). For the superposed magnetic flux from star and disk we 
have 



*(r, z) = ^ ,d \ (\/r 2 + (z d + z) 2 -{z d + z 

+ *0, 



(r 2 + (z d + z 



2^3/2' 



(7) 



with the stellar and disk magnetic flux \I/o,* and ^o,d- The poloidal magnetic field follows 
from the derivatives B r = —(l/r)(d^>/(dz) and B z = (l/r)(d^f /(dr), properly calulated in 
the staggered mesh in order to obtain a numerically divergence- free and force-free initial field 
structure. The dimensionless disk thickness z d with (z d + z) > for z < is introduced in 
order to avoid kinks in the field distribution. Several combinations of both field components 
were investigated, parameterized by combinations of ^/ ,d and \l/ ,*- Figure [2] shows three 
examples for the initial field configuration for different strength and alignment of the fiel d 
components. Sim ilar field geometries have been discussed already by lUchida &: Low! (1198lh . 
Matt et al.l (120021 ) did consider similar configurations, however, superimposing a stellar dipole 
with a weak vertical disk field. 

In order to allow for a clear comparison between all our runs, the respective magnetic 
field component profiles and wind density component profiles are identical. 



4. Results and discussion 

We now discuss a number of MHD jet formation simulations covering a wide parameter 
range (see Tab. [1]) concerning both the magnitude of the disk and stellar wind mass load 
and magnetic flux, respectively. The results presented here are preliminary in the sense that 
not all simulation runs could be performed over time scales sufficiently long enough for the 
MHD flow as a whole to reach the grid boundaries or to establish a stationary state. This is 
due to the comparatively large physical grid size in combination with steep gradients in the 
disk wind parameters which in general allow only for a weak outflow from large disk radii. 
In particular the steep decline of the stellar field in combination with a reasonable disk mass 
loss rate is numerically problematic (see below). Therefore, for a true comparison between 
different runs it is essential to also take into account the dynamical state of the outflow. 

In order to check for resolution issues, we have also run a set of simulations with twice 
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Fig. 3. — Time evolution of simulation run A4a. Poloidal magnetic field distribution at time 
steps t = 0, 50, 200, 500, 1000, 2000, 2500, 3000, 3600 (from top left to bottom right). Density 
gray scale as indicated. Poloidal field lines as contour levels of the magnetic flux, — \P(r, z) = 
3.0 ,2.0 ,1.6 ,1.3, 1.0, 0.8, 0.6, 0.4, 0.2, 0.12, 0.064, 0.032, 0.016, 0.008, 0.004, 0.002, 0.001, 
0.0005, 0.0002, 0.0001, (dashed), *(r, z) = 0.0002, 0.0005, 0.001, 0.002, 0.004, 0.008, 0.016, 
0.032, 0.064, 0.12, 0.2, 0.3, 0.4, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0 (solid). 
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Fig. 4. — Time evolution of the density at z — 44, r = 44 for simulation run A4a. 




Fig. 5. — Radial profiles of density, magnetic flux, axial magnetic field, axial velocity, and 
radial velocity (from top to bottom) at z = 55 for simulation run A4a at t — 3000. Note 
the neutral field line (magnetic field reversal) at r ~ 33 and the magnetic flux reversal at 
r ~ 53. 
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Fig. 6. — Poloidal magnetic field evolution during two flares around t = 176 and t = 210. 
Solid and dashed lines indicate the direction of total magnetic flux of the superposed dipolar 
and disk magnetic field components. Shown are time steps: 1700, 1760, 1790, 1810, 1840, 
1890, 2080, 2100, 2130, 2140, 2180, 2250 (from top left to bottom right) of simulation A4a. 
Contour levels as in Fig. [31 
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the resolution, without seeing significant differences (runs A4b, A5, A6). Thus, since we 
are interested in many parameter runs evolving for a long time, we mainly concentrate on 
the lower resolution simu lations. We describe first the general evolution of jet format ion in 
our simulations (see also fbuved fc Pudritzlll997l : [Fendt fc Cemeliidbooi iFendtlhood n and 
later discuss specific features. 



4.1. Overall evolution of an outflow 

Figure [3] shows how the field structure evolves in time for the simulation run A4a. In 
this case the stellar dipolar field and the disk field are aligned. Thus, the X-point is initially 
along the rotational axis. We have stopped this simulation after 3600 Keplerian rotations 
at the inner disk radius. This corresponds to 6 rotations at the outer disk radius and makes 
it clear that the outer parts of the disk wind have not yet completely evolved into a quasi- 
stationary state. This is a generic problem of all disk-jet simulations published so far as soon 
as large disk radii are considered. 

The evolution during the very first time steps from the initial shows that above the 
emerging outflow the initial steady state corona is still present. The X-point which was 
initially located at z < 10 along the rotational axis, moves upwards to z ~ 20 (for t = 50) 
and z ~ 40 (for t = 200) and is later swept out of the computational domain for t > 600 
(see top of Fig. [3] and discussion end of Sect. 4. 2). 

Since the initial condition is still kept in steady state, artificial dynamical re-configuration 
is prevented. This demonstrates that the numerical resolution is sufficient for our problem 
also in the out parts of the grid. 

The next dominant feature is observed still at early stages of about ten rotations. The 
initial dipolar field breaks up due to the magnetic pressure of the toroidal magnetic field 
induced by differential rotation between star and disk. The magnetic pressure gradient drives 
this outflow. Furthermore, an intermediate axial jet is launched due to the re-arrangement 
of the initially hydrostatic corona to a new dynamical equilibrium state. In fact, as the 
magnetic field along the axis is squeezed by lateral dynamical pressure, the material is 
accelerated along the rotational axis. After the dipole is broken up, there is few direct 
magnetic connection between star and disk anymore and the differential rotation induced 
toroidal field decreases. 



The outflows from star and disk continue to grow gaining higher kinetic energy and 
momentum. MHD self-induction of toroidal magnetic field leads to collimation and magnetic 
acceleration. Compared to pure disk winds (jOuyed fc Pudritzlll997l ; iFendt fc Cemeljidl2002l ; 



Fendtll2006l ) the outflow is clearly less collimated as being de-collimated by the central stellar 
wind. 



At intermediate time scales quasi- stationary states may emerge. This is demonstrated 
for example in simulation A4a in Fig. H] by showing the time evolution of the density at point 
(z = 44, r = 44). Two plateaus are clearly seen at different density levels indicating two 
quasi- stationary states of the flow evolution. One is from t ~ 1000 till t ~ 1500, the other 
from t ~ 2300 — 3000 when the outflow is slowly re-adjusting from the flaring events (see 
below). Note, however, that at this time the outer disk has rotated only about a fifth of an 
orbit. Thus, the field and flow above the outer disk will further evolve in time and again 
disturb the enclosed structure in quasi-steady state. We observe that over even longer time 
scale such quasi-stationary states may be reached (and be disturbed) again and probably 
again and again. We believe that this feature is due to the still ongoing evolution of the 
outer disk wind. The final states of some representative simulation runs are shown in Fig. [TBI 
and [T41 We will discuss them in Sect l4.4.3l below. 

As seen already from the initial field configuration, the simulations of the anti-aligned 
field configuration (as e.g. A4a) reveal a change of sign in the magnetic flux distribution 
from star and disk (negative magnetic flux is indicated by dashed contours in Fig. [5]). The 
location of flux reversal is accompanied with a concentration of contour lines in our figures 
which may be confused with the existence of a shock. This is however an artifact of the 
choice of contour levels for the magnetic flux. In fact, it makes sense to follow the same 
field lines emerging from the stellar surface also back into the disk surface. Since the strong 
gradient of the stellar dipolar field, magnetic flux levels aroung \1/ = are concentrated. 
Figure [5] shows that the radial profiles of density, velocity and poloidal field accross the flux 
reversal. Note that the location of magnetic field reversal (the neutral field line) is at a 
smaller radius r = 33 compared to the magnetic flux reversal at about r = 53. Within the 
magnetic field reversal magnetic flux is accumulated. Beyond the field reversal the magnetic 
flux decreases and then becomes negative. 



4.2. Reconnection and large scale flares 



Our version of ZEUS code has physical magnetic diffusivity implemented (see lFendt fc Cemeljic 



(120021 ) for explanation and tests). This allows to consider reconnection processes. 

We observe reconnection flares along some of the outflows (see Fig. [6]). These flares are 
similar to coronal mass ejection. They rapidly evolve and propagate along the neutral field 
line. Once formed, reconnection islands (or rather "tori" in our axisymmetric setup) propa- 
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gate across the jet magnetosphere within a few rotation times and leave the computational 
domain. The flares typically expand and reconnect within 70 time units equivalent to 70 
orbital periods of the inner disk (and 70 stellar rotations). We may also naively measure the 
flare propagation speed by the proper motion of the field lines observed in the simulation. 
This gives an average "flare propagation speed" of about unity, as they travel 80 spatial 
units in 70 time units - the Keplerian speed at their foot-point. The flares propagate along 
the neutral field line, and leave the physical grid at a radius corresponding to about 7-10 
AU from the axis. 

In the following we mainly concentrate on simulation A4a. Figure [11] is a magnification 
of the inner area of Figj6] and demonstrates that the flares are actually from a Y-point above 
the disk, located at (z ~ 2, r ~ 4) Note that the stellar magnetosphere remains closed also 
for foot-points along the disk with r < 3 (at t ~ 1700), respectively r < 2 (t ~ 1700), beyond 
the co-rotation radius. In this case the flare eruptions are launched close to the inner disk 
area, but not at the inner disk radius. 



Goodson et al.l (119991 ) also observed flares or reconnection events in their simulations of 
a dipolar magnetosphere connected top a surrounding accretion disk. Since the evolution of 
the disk structure was treated in their simulation, they were able to observe time-dependent 
accretion along the equatorial plane. The frequent ejection events of AU-sized knots along the 
rotational axis were correlated with a time-variation of the accretion rate. This is different 
in our case, where the reconnection/flares seem to be triggered by the evolution of the 
outer disk wind. As mentioned, during the run-time of our simulations (up to 3600 inner 
disk rotations) the corona (or disk wind respectively) above the outer disk has not yet 




Fig. 7. — Integrated mass flux in axial direction across the upper ^-boundary versus time. 
Note the change of mass flux of about 10 — 50% during the flare events. The high mass flux 
for t < 500 indicates sweeping off of the initial hydrostatic corona. 
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Fig. 8. — Profile of the outflow momentum along the z-boundary at time t = 
1400, 1600, 1800, 1900, 2100, 2200, 2500, 3000. Note the re-distribution of the main mass flow 
channel from radius r = 60 to r = 40. 
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Fig. 9. — Profile of the axial outflow velocity along the z-boundary at time t = 
1400, 1800, 1850, 2200, 2300, 3000. 
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reached a steady dynamical statqj. Thus, in our case, as the outer disk outflow evolves, the 
cross-jet pressure equilibrium is changed and forces the inner magnetic field configuration to 
adjust accordingly. Since Goodson et al did not specify the distribution and magnitude of 
the magnetic diffusivity in their simulations, it is difficult to compare the results in detail. 
Another difference is the duration of the simulation. While our physical grid size is only 
marginally larger, the time scale of our simulations is substantially longer (3600 inner disk 
rotations corresponding to about 36,000 days compared to 130 days). This is particularly 
important when considering the dynamical evolution above the outer part of the disk and 
the long-term behavior. 

The flare events observed in our simulations are accompanied by a temporal change in 
outflow mass flux, velocity, or momentum, respectively. Figure [7] shows the mass loss rate in 
axial direction integrated across the upper r-boundary versus time. During the first flare we 
see a 10%-variation in the mass flux followed by a sudden decrease of mass flux by a factor 
of two during flare two. Figures [8] and [9] show the profiles of jet momentum and poloidal 
velocity across the upper r-boundary. We see that during flaring these profile indicate a 
re-arrangement of momentum and velocity profile across the jet. Before the flare, the radial 
jet momentum profile is broad (see Fig. [H] upper left for t = 1400). In fact, the profile 
(Fig. El t = 1400) remains very similar for several 100 rotations before until the flare event 
starts. After the flare has passed the grid, the jet momentum profile is concentrated within 
a cylindrical sheet of radius r ~ 35 and thickness Ar ~ 5 (Fig. [8] lower right at t — 2300). 
Similar to the situation before the flare, the jet momentum profile for time steps after the 
flare (t > 2300) look almost identical (see bottom sub-figures for t = 2600 and t = 3000). 

This behavior is somewhat mirrored in the poloidal velocity profile. What is interesting 
for shock formation in the asymptotic jet, is a temporal change in jet velocity at certain 
jet radii (see below). The maximum velocities in the outflow reach about two times the 
Keplerian speed at the inner disk radius and vary by a factor of two. Along the neutral field 
line the outflow velocity is fast-magnetosonic (see Figs. [11] and [14] where the fast surface is 
indicated). 

As an estimate for the reconnection time scale we may apply the Sweet-Parker approach 
with tsp ~ yj ta Tdifr with the dynamical (Alfven) time scale ta = l/v\ and the diffusive time 
scale T^ff = l 2 /v- With our model of turbulent magnetic diffusivity r] = r^op 1 ^ 3 we find 
r SP ~ / 3 / 2 (47r) 1 / 4 r] " 1 ^ 2 p 1 / 12 . For r] = 0.01, and p ~ 0.02 at the reconnection area with I ~ 2 
(for simulation A16) that tsp ~ 30. This is similar to the duration of the reconnection flare 
we observe between t = 600 and t = 650. 



3 However, the same argument holds for the Goodson et al. work 
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Note an interesting feature in both of the general model setups (aligned and anti-aligned 
field geometry) concerning the reconnection geometry. Figures [11] and [12] show the inner 
structure of the star-disk magnetosphere. These simulations (A04a, A16) were launched 
from differently aligned field geometries (aligned/anti-aligned disk-stellar field). In the anti- 
aligned case the initial X-point is along the axis. However, as the magnetic field evolves, 
the dipole expands and sweeps off the initial field along the axis. A new X-point evolves by 
disruption of some of the closed dipolar field loops. The new X-point is located above the 
disk close to the inner disk radius (see above). This situation is not very different from the 
the aligned case for which the initial X-point, located along the equatorial plane, moves up to 
a certain (small) height above the disk. In summary, our simulations show that no magnetic 
X-point remains along the equatorial plane as e.g. assumed in the Shu et al. X-wind model. 
Instead we find from both initial co nfigurations an X-point above the disk (sometimes also 
called Y-point, (IFerreira et al.ll2006l )). located at time-averaged distances (z ~ 4, r ~ 2) or 
(z ~ 8,r ~ 2), respectively. 



4.3. Flaring events as jet knot generator? 

The generation of knots in protostellar jets is a long-standing puzzle. In particular, it is 
unclear whether the shocks/knots are triggered by the internal engine or due to interaction 
with the ambient medium. A strong arg ument for the first po ssibility is the existence of 



some perfectly symmetric jets as HH212 (IZinnecker et al.lll998l ). however, the majority of 



jet sources looks asymmetric. Internal shocks along the jet flow may be cau sed e.g. by 



a tim e-dependent velocity variation of the material injected into the jet (e.g. iRaga et al. 



( 120071 )). Even if we have observational support for the idea of the jet knots being triggered 
by the central engine, we do not know how the central engine does that. The time scale 
derived from typical knot separations d^ QOt and velocities fknot correspond to 



iknot = 10 yrs 



^knot \ f d 



300km/ S y V100AU 



^knot 



which is clearly different from the Keplerian time scale close to the inner disk edge from 
where the jets are launched. 

Considering the ejection of large-scale flares and the follow-up re-arrangement of outflow 
density and velocity distribution in our simulations, one is tempted to hypothesize that the 
origin of knots is triggered by such flaring events. 

In our simulations (A4a, A16) the time scale of flare generation is about 500-1000 
rotational periods of the inner disk, corresponding to about 30-60 years (assuming an inner 
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disk rotation period of, say two times the co-rotation period of a typical TTauri star of 
about 10 days). The variation in the hydrodynamic parameters lasts for about 30-40 inner 
disk rotation times (respectively about < 400 days for a 10 day stellar rotation period). 
The further evolution and generation of further flaring events is of course not known as it 
is beyond our computation time. The essential point, however, is that we detect a long 
time scale, substantially longer than the typical dynamical time scales at the jet formation 
area. That time scale is surprisingly similar to the time scale defined by the observed knot 
separation and velocity. Of course, it is to early to draw firm conclusions from such a tentative 
agreement. The reconnection time-scale is governed by the magnetic diffu sivity for which we 



have applied our self-consistent model of turbulent magnetic diffusivity (IFendt fc Cemeljic 



20021 ). 



It is left to further simulations of the asymptotic collimated jet beam to check whether 
the detected variation in outflow speed and mass flux (Figs. [71IH1I91 is sufficient in order to 
generate strong internal shocks comparable to the observed knots. 



4.4. Collimation degree, mass loss rate and field alignment 

In this section we co mpare the overall collimation behavio r of the different simulation 



runs. In previous studies (IFendt fc Cemeljid 120021 ; iFendtl 120061 ) we quantified the degree of 



outflow collimation ( by comparing the mass flow rates in axial and lateral direction, 

( = Ml= 2n rpVzdr (9) 
M r 27rr max pv r dz. 

In Tab. [T] we provide mass flow rate within four differently sized volumes, M Zj j, M r ,*i consider- 
ing different sub-grids of the whole computational domain, i.e. cylinders of radius and height 
(fmax x -2max) = (12 x 12), (24 x 24), (43 x 43), and (76 x 76). We calculate the mass fluxes 
ratio normalized to the size of the area threaded by the mass flows in r- and z-direction, 
Cij C2, C35 Ci- The ratio of grid extension (r max / 'z max ) as displayed in our figures converts into 
a ratio of cylinder surface areas of A Z /A T ~ 0.5(r max /z max ). As a general measure of collima- 
tion degree, we give an average value < ( > which may consider also the dynamical state of 
the simulation run. 

Compared to our previous work on pure disk wind collimation, where the relative mass 
flux of the collimating disk wind is a sufficient quantitative measure of collimation the sit- 
uation in the present setup is not exactly the same. Naturally, the stellar wind mass flux 
emerges close to the outflow axis and will tend to stay close to it. Thus, a high stellar mass 
flux will naturally cause a more collimated mass flow. 
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In order to compare the collimation degree within the large volume, it is essential to 
check the evolutionary state of the flow. In all of our simulation runs the initial corona has 
completely swept out of the grid. In this case, all of the mass flux we measure has been 
launched by the disk wind. Exceptions are simulations A10, A13, and A12, where a relict 
of the initial bow shock are still visible in the outer layers. The mass loss rates in these 
examples have to taken with care, at least the values from the larger volumes. Run A16 is 
evolving very slowly in the outer part and it was not possible within reasonable CPU time 
to evolve the disk wind to larger distances from the disk surface. 



4-4- 1- Collimation along the outflow 

The degree of collimation derived from the simulations is in general different for the the 
different volumes, Cij C2, C35 Ca- 
in most cases, the degree of mass flux collimation increases along the flow which is ex- 
actly the signature of MHD self-collimation. The outflow starts as un-collimated disk/stellar 
wind, and reaches the outer grid boundaries as collimated beam. 

A good example is simulation A7, where the collimation of mass flow changes from 
ration ~ 0.3 for the inner region to ~ 7 in the outer parts. The same arguments hold for 
simulations A8, A15, A14. A similar behavior is seen e.g. in simulations A10, A13, A14, 
however, the collimation in mass flux around the largest volume is low. In these cases the 
outermost flow structure has not yet evolved into a steady state and either parts of the initial 
bow shock or the initial corona are still in the computational domain. 

A counter example is simulation A12, which stays un-collimated in mass flux, although 
the magnetic field distribution looks collimated. Simulation A12 is exceptional for its 
high mass flux launched from the disk surface. Thus, the disk wind starts super-Alfvenic 
and super-seed the magnetosonic speed quickly. The standard Blandford-Payne magneto- 
centrifugal acceleration is not very efficient in this case, which also results in only little in- 
duction of toroidal magnetic field and, thus, Lorentz force, see FigJTOl The radial component 
of the perpendicular Lorentz force component is negative, i.e. directed inwards, however, 
too weak to balance the strong inner centrifugal force, resulting in a weak collimation of the 
outer part. 

Table [1] shows also the time evolution of the collimation degree. Example A15 demon- 
strates how the collimation degree progresses in time. From earlier times (t=550) to later 
times (t = 2100) collimation increases as the outflow evolves into a new dynamical state over 
the whole numerical grid. The same hold for example A14 and others. Examples A7, A8, 
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A9 are calculated for a high stellar wind mass load, thus the time-evolution of the outer disk 
wind does not play a big role concerning collimation. This outflows reach a high degree of 
collimation early in time. 

4-4-2. Collimation and mass loss rates 

The mass loss rates from star and disk are prescribed as boundary condition. A higher 
mass load in the stellar wind component would naturally result in a higher degree of colli- 
mation in the mass flux as this mass is just launched closer to the outflow axis. 

Simulations A2 and A3 are good examples. The stellar wind mass loss is 2-3 times the 
disk wind mass flux. The resulting mass flux ratio is similar, thus this outflow is rather well 
collimated hydrodynamically. However, in spite the quite good collimation in mass flux, the 
magnetic field structure is not collimated with almost (spherical) radially expanding poloidal 
field lines (see Fig. [TS]). The field structure (and thus the poloidal velocity field) has a conical 
shape and is thus un-collimated. Similar as for simulation A12 (see above) the disk wind is 
too weak in order to collimate the stellar wind to a high degree. 

Simulations A12 and A13 have the same star-disk magnetic field profile, however with 
a different mass load by a factor of ten. Run A12 reaches the same dynamical evolutionary 
state as A13, but earlier at t — 300 instead of t = 620. The less magnetized outflow A13 is 
stronger collimated. Run A14 has the same mass load as A13, but the double magnetic flux 
resulting in similar collimation degree. 

Simulations A10 and A13 differ only by the stellar wind mass load (factor two), however 
this mass load is low and does not result in a variation of the overall degree of collimation. 
Note that for both simulations, the outer area of the outflow has not yet evolved in a steady 




Fig. 10. — Perpendicular Lorentz force (top) and centrifugal force (bottom) across simulation 
A12 at height z = 50 at time t = 380. 
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state and a relict from the initial bow shock are still visible. 



4-4-3- Collimation and magnetic field alignment 

From the technical point, the simulation of the case where magnetic dipole and disk 
magnetic field are anti-aligned (A3, A4, A10, A13, A12, A14) are computationally less 
expensive and were therefore able to evolve to considerably longer time. 

For the aligned cases (A9, A8, A7, A16, A15, A2) only the less magnetized evolve 
into a quasi stationary state in the outer (disk) outflow in reasonable time. We mention in 
particular run A15 with a reasonably strong disk magnetic flux which enables a collimating 
disk wind (compare to A16 and A7). 

Simulations A9, A8 and A7 demonstrate that a low mass flux disk wind mass flux 
cannot evolve into a collimated disk wind even if the magnetic flux is relatively high if it is 
dominated by a strong central stellar jet. 

The most promising model setup in order to explain strong stellar jets from a star-disk 
magnetosphere are those with relatively strong disk wind and disk magnetic flux. The stellar 
wind dominated simulations may give a high degree of collimation, however they collimate 
to too small radii. Stellar magnetic flux dominated simulations tend to stay un-collimated. 

5. Summary 

We have performed axisymmetric MHD simulations of jet formation from accretion 
disks surrounding a magnetised star. The simulations start from an equilibrium state of the 
star-disk corona (hydrostatic density distribution, force-free field). Our physical grid size is 
(80 x 80) inner disk radii corresponding to about (14 x 14) AU for r; ~ 10 R Q . Disk and 
stellar surface are taken as a time-independent boundary condition for the outflow mass loss 
rates and the magnetic flux profile. 

The major goal of this paper was to investigate the long-term interrelation between a 
stellar dipolar field and a disk field and how that affects the outflow collimation. Certain 
combinations of stellar versus disk magnetic flux and field directions were considered namely 
the cases of an aligned resp. anti-aligned magnetic field direction. 

The stellar magnetic field has an important impact on the jet formation process by 
providing additional magnetic flux, additional central (magnetic) pressure component, excess 
angular momentum in the jet launching region. It might disturb the outflow axisymmetry, 
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Fig. 11. — Poloidal magnetic field distribution (thin lines) during flare evolution (run A04). 
Contour levels as in Fig. Simulation time steps t=2000, 2060, 2070, 2130, 2140, 3000. 
(from top left to bottom right). Alfven and fast surface are shown by thick lines. Note the 
super- Alfvenic /super-fast stellar wind. 
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Fig. 12. — Poloidal magnetic field distribution (thin lines) during flare evolution (run A16). 
Contour levels as in Fig. O Simulation time steps t= 600, 620, 630, 640, 670, 760. (from 
top left to bottom right). Alfven and fast surface are shown by thick lines. Note the super- 
Alfvenic /super-fast stellar wind. 
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Fig. 13. — Poloidal magnetic field lines at the time when the simulation was stopped. Sim- 
ulation A9 (t=740); A8 (t=1000); A7 (t=1400); A16 (t=1350) A15 (t=2100); A2 (t=700) 
(from top left to bottom right). Thick and medium thick lines indicate the location of the 
Alfven and the magnetosonic surface. Contour levels and grey scale as in Fig. [3j 




Fig. 14. — Poloidal magnetic field distribution at the time when the simulation was stopped. 
Simulation A3 (t=500); A4a (t=3200); A10 (t=550); A13 (t=620); A12 (t=380); A14 
(t=660) (from top left to bottom right). Contour levels and grey scale as in Fig. [3J 
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Table 1: Summary of simulation runs ID. Shown are disk and stellar mass loss rates, M*, 
corresponding magnetic fluxes, ^o,d, ^o,*, and the physical time step r, when the collimation 
degree has been calculated. Mass loss rates M zi , M ri in z- and r-direction are integrated along 
three sub-grids (r ijmax x z i:iaax ) = (11.6 x 11.6), (23.6 x 23.6), (43.4 x 43.4), (76.2 x 76.2), (for 
the (80.0 x 80.0) grid). The average degree of collimation <C> is defined by the relative 
mass fluxes in z and r-direction, normalized by the area threaded. Simulations labeled with 
a '*' were run on a smaller, higher resolution (40.0 x 40.0) grid. 
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and/or trigger a time- variation in outflow rate in the case of an inclined central stellar 
magnetosphere. 

Our results of MHD simulations of a superposed stellar & disk magnetosphere in gen- 
eral demonstrate the de-collimation of the disk jet by the central stellar magnetosphere, 
respectively the collimation of the stellar wind by the surrounding disk jet. 

The interplay between disk wind and stellar wind, may result in flares like coronal mass 
ejections, generating a variation in the mass loading and velocity of the asymptotic jet. We 
find variations in mass load by a factor of four and in velocity by a factor of two. The 
time scale for such flaring events is long, i.e. several hundred inner disk rotations and of 
the order of 10.000 days. We therefore hypothesized whether such flaring events and the 
corresponding hydrodynamic variations may be responsible for generating internal shocks in 
the asymptotic jet visible as knots. 

In summary we conclude that strong protostellar jets must be generated by collimating 
disk winds with reasonable mass load. If the two-component system is dominated by the 
stellar wind, collimation is either too weak (for low mass load) or too high (for high mass 
load). It is only a disk jet which can provide a collimated outflow over a reasonable range 
in radius. 

I thank the LCA team and M. Norman and D. Clarke for the possibility to use the 
ZEUS code. Comments and suggestions by an unknown referee which helped to improve the 
presentation of this work are acknowledged. 
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